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Here  considered  is  the  mathematical  analysis  and  numerical  computation  of  the  electro¬ 
magnetic  wave  scattering  by  multiple  cavities  embedded  in  an  infinite  ground  plane. 
Above  the  ground  plane  the  space  is  filled  with  a  homogeneous  medium,  while  the  interi¬ 
ors  of  the  cavities  are  filled  with  inhomogeneous  media  characterized  by  variable  permit¬ 
tivities.  By  introducing  a  new  transparent  boundary  condition  on  the  cavity  apertures,  the 
multiple  cavity  scattering  problem  is  reduced  to  a  boundary  value  problem  of  the  two- 
dimensional  Helmholtz  equation  imposed  in  the  separated  interior  domains  of  the  cavities. 
The  existence  and  uniqueness  of  the  weak  solution  for  the  model  problem  is  achieved  via  a 
variational  approach.  A  block  Gauss-Seidel  iterative  method  is  introduced  to  solve  the  cou¬ 
pled  system  of  the  multiple  cavity  scattering  problem,  where  only  a  single  cavity  scattering 
problem  is  required  to  be  solved  at  each  iteration.  Numerical  examples  demonstrate  the 
efficiency  and  accuracy  of  the  proposed  method. 
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1.  Introduction 

The  phenomenon  of  electromagnetic  scattering  by  cavity-backed  apertures  has  received  much  attention  by  both  the  engi¬ 
neering  and  mathematical  communities  for  its  important  applications.  For  instance,  the  radar  cross  section  is  a  measure  of 
the  detectability  of  a  target  by  a  radar  system.  Deliberate  control  in  the  form  of  enhancement  or  reduction  of  the  radar  cross 
section  of  a  target  is  of  no  less  importance  than  many  radar  applications.  The  cavity  radar  cross  section  caused  by  jet  engine 
inlet  ducts  or  cavity-backed  antennas  can  dominate  the  total  radar  cross  section.  A  thorough  understanding  of  the  electro¬ 
magnetic  scattering  characteristic  of  a  target,  particularly  a  cavity,  is  necessary  for  successful  implementation  of  any  desired 
control  of  its  radar  cross  section,  and  is  of  high  interest  to  the  scientific  and  engineering  community.  Another  example  is  that 
the  cavity  can  be  used  to  model  cracks  or  holes  in  metallic  surfaces  such  as  aircraft  wings.  These  cracks  or  holes  would  be 
invisible  to  a  visual  inspection  but  may  be  revealed  by  understanding  the  scattering  characteristics  of  the  cavity.  As  an  in¬ 
verse  problem,  the  mathematical  model  can  serve  as  a  predictor  of  the  scattering  of  electromagnetic  waves  by  the  cavity  for 
use  in  non-destructive  testing.  Besides,  this  work  is  also  motivated  by  the  study  of  the  optimal  design  problems  of  the  cavity, 
where  one  wishes  to  design  a  cavity  structure  that  reduces  or  enhances  the  radar  cross  section. 

For  the  time-harmonic  analysis  of  cavity-backed  apertures  with  penetrable  material  filling  the  cavity  Interior,  we  men¬ 
tion,  for  example,  works  by  Jin  et  al.  [23,24,27],  Wood  and  Wood  [37],  and  references  cited  therein.  Mathematical  analysis 
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of  the  cavity  scattering  problem  including  the  scattering  from  overfilled  cavities,  where  the  cavity  aperture  is  not  planar  and 
may  protrude  the  ground  plane,  can  be  found  in  Ammari  et  al.  [1-3],  Bao  et  al.  [9],  Li  et  al.  [26],  Van  and  Wood  [28-33],  Wood 
[36].  Mode  matching  based  analytical  approaches  are  developed  in  Bao  and  Zhang  [13],  and  Bao  et  al.  [10]  for  solving  the 
electromagnetic  scattering  problem  involving  large  cavities.  Much  research  has  been  devoted  to  solving  tbe  cavity  scattering 
problem  by  various  numerical  methods,  including  finite  element,  finite  difference,  boundary  element,  and  hybrid  methods. 
See,  for  example,  Bao  and  Sun  [11],  Du  [15,16],  Huang  and  Wood  [19],  Huang  et  al.  [20],  Wang  et  al.  [34,35],  Zhao  et  al.  [39] 
and  Zhang  et  al.  [40[.  To  the  best  of  our  knowledge,  all  the  known  results  in  the  open  literature  follow  the  model  of  a  single 
cavity,  which  clearly  limits  the  practical  application  of  the  model  problem  in  industiy  and  military.  A  major  challenge  in  a 
multiple  cavity  model  is  how  to  fully  capture  the  Interactions  among  separate  cavities.  This  paper  aims  to  extend  the  single 
cavity  model  to  the  more  general  multiple  cavity  model,  analyze  and  develop  numerical  methods  for  the  associated  scatter¬ 
ing  problem.  Some  preliminary  numerical  results  are  announced  in  Li  and  Wood  [25]. 

Our  approach  calls  for  the  development  of  a  boundary  condition  over  the  cavity  apertures  based  on  Fourier  transform.  The 
boundary  condition  is  nonlocal  and  transparent  which  connects  the  fields  in  all  individual  cavities.  By  using  this  boundary 
condition,  we  reduce  the  multiple  cavity  scattering  problem  into  a  boundary  value  problem  of  the  two-dimensional  Helm¬ 
holtz  equation  imposed  in  the  interiors  of  the  cavities.  The  existence  and  uniqueness  of  the  weak  solution  of  the  associated 
variational  formulation  for  the  model  problem  Is  achieved.  A  block  Gauss-Seldel  Iterative  method  Is  Introduced  to  solve  the 
coupled  system,  where  only  a  single  cavity  scattering  problem  is  required  to  be  solved  at  each  iteration,  rendering  all  effi¬ 
cient  single  cavity  solvers  applicable  to  the  multiple  cavity  scattering  problem.  Two  numerical  examples  are  presented  to 
show  the  efficiency  and  accuracy  of  the  proposed  method.  Numerical  methods  for  multiple  obstacle  scattering  problems 
can  be  found  In  Grote  and  Kirsch  [17],  Huang  et  al.  [18],  Jiang  and  Zheng  [22],  and  references  therein.  Recently,  Bonnetier 
and  Triki  [14]  and  Babadjian  et  al.  [8]  studied  the  enhancement  of  electromagnetic  fields  by  interacting  subwavelength  cav¬ 
ities.  More  general  results  could  be  found  in  Ammari  et  al.  [6]  on  the  existence  of  resonances  by  using  layer  potential  tech¬ 
niques.  We  also  refer  to  Ammari  et  al.  [4,5,7]  for  related  work  on  developing  imaging  functions  to  reconstruct  perfectly 
conducting  cracks  and  holes. 

We  organize  the  paper  as  follows.  Section  2  presents  a  mathematical  model  for  the  single  cavity  scattering  problem;  the 
variational  formulation  is  presented  for  an  equivalent  boundary  value  problem  by  using  the  transparent  boundary  condition; 
the  uniqueness  and  existence  of  the  solution  are  examined.  Section  3  deals  with  two  cavities.  The  major  new  ingredient  is 
the  introduction  of  a  novel  transparent  boundaiy  condition  which  reduces  the  double  cavity  problem  to  two  coupled  single 
cavity  problems.  The  results  are  extended  to  the  general  multiple  cavity  scattering  problem  In  Section  4.  In  Section  5,  we 
examine  issues  related  to  the  numerical  implementation.  Experiments  on  two  example  scenarios  are  performed  and  shown 
to  be  efficient  and  accurate.  The  paper  is  concluded  in  Section  6. 

2.  One  cavity  scattering 

In  this  section,  we  study  a  mathematical  model  for  a  single  cavity  scattering  problem,  and  discuss  the  existence  and 
uniqueness  of  the  solution  based  on  its  variational  formulation.  This  section  is  intended  to  introduce  the  background  for 
the  cavity  scattering  problem  and  serve  as  a  basis  for  the  two  cavity  and  the  general  multiple  cavity  scattering  problems. 

2.1.  A  model  problem 

We  focus  on  a  two-dimensional  geometry  by  assuming  that  the  medium  and  material  are  invariant  in  the  z-direction. 
Throughout,  the  medium  is  assumed  to  be  non-magnetic  and  has  a  constant  magnetic  permeability,  i.e.,  ft  =  fig,  where  fig 
is  the  magnetic  permeability  of  vacuum.  The  electromagnetic  property  of  the  medium  is  characterized  by  the  dielectric  per¬ 
mittivity  e. 

As  shown  in  Fig.  1,  an  open  cavity  fi,  enclosed  by  the  aperture  F  and  the  wall  S,  is  placed  on  a  perfectly  conducting  ground 
plane  F"^.  Above  the  flat  surface  {y  =  0}  =  F  u  F*^,  the  medium  is  assumed  to  be  homogeneous  with  a  positive  dielectric  per¬ 
mittivity  fio.  The  medium  inside  the  cavity  Q  is  inhomogeneous  with  a  variable  relative  dielectric  permittivity  e{x,y).  Assume 
further  that  e{x,y)  e  L°°(Q),  Res  >  0.  Ime  >  0. 

Due  to  the  uniformity  in  the  z-axis,  the  scattering  problem  can  be  decomposed  Into  two  fundamental  polarizations;  trans¬ 
verse  magnetic  (TM)  and  transverse  electric  (TE).  Its  solution  then  can  be  expressed  as  a  linear  combination  of  the  solutions 
to  TM  and  TE  problems.  In  this  work,  we  focus  on  the  TM  polarization.  The  method  can  be  extended  to  the  TE  polarization 
with  obvious  modifications. 


IX  ,,  IX 


n 


s 

Fig.  1.  The  problem  geometry  for  a  single  cavity  scattering  problem.  An  open  cavity  Q,  enclosed  by  the  aperture  F  and  the  wall  S,  is  placed  on  a  perfectly 
conducting  ground  plane  F*^. 
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For  the  TM  polarization,  the  magnetic  field  is  transverse  to  the  invariant  direction.  The  incident  and  the  total  electric  fields 
are  parallel  to  the  invariant  dimension.  By  the  perfectly  electrical  conductor  condition,  the  total  field  u  vanishes  on  the  cavity 
wall  S  and  the  ground  plane  F'  except  over  the  aperture  F.  The  time-harmonic  Maxwell  equations  are  reduced  to  the  two- 


dimensional  Helmholtz  equation 

Au  +  k^u  =  0,  in  flu  R^, 

(2.1) 

together  with  the  homogeneous  Dirichlet  boundary  condition 

u  =  0,  on  F*^  U  S. 

(2.2) 

Here  =  co^efig,  where  to  is  the  angular  frequency  and  k  is  known  as  the  wavenumber. 

Let  an  incoming  plane  wave  u'  =  ei^colxsinfl-ycoso)  incident  on  the  cavity  from  above,  where  9  e  (-71/2,  n/l)  is  the  angle  of 
incidence  with  respect  to  the  positive  y-axis,  and  Kq  =  oj^Sofig  is  the  wavenumber  of  the  free  space. 

Denote  the  reference  field  as  the  solution  of  the  homogeneous  Helmholtz  equation  in  the  upper  half  space: 

=  0,  in  IR2  ,  (2.3) 

together  with  the  boundary  condition 

=  0  on  F"  U  F.  (2.4) 

It  can  be  shown  from  (2.3)  and  (2.4)  that  the  reference  field  consists  of  the  incident  field  u'  and  the  reflected  field  uF 

yref  ^  iji  + 

where  u""  =  _eiKo(xsino+ycos0)_ 

The  total  field  u  is  composed  of  the  reference  field  and  the  scattered  field  uF 
u  =  +  u* . 

It  can  be  verified  from  (2.1)  and  (2.3)  that  the  scattered  field  satisfies 


Au^  +  KqU*  =  0,  in  [R^. 

In  addition,  the  scattered  field  is  required  to  satisfy  the  radiation  condition 

fdu^ 


limyp 


dp 


-ijcou*  =0,  /i  =  |(x,y)|. 


(2.5) 


(2.6) 


To  describe  the  boundary  value  problem  and  derive  its  variational  formulation,  we  need  to  introduce  some  functional 
space  notation.  For  u  e  L^(F'^  u  F),  which  is  identified  with  L^(R),  we  denote  by  Ci  the  Fourier  transform  of  u  defined  as 


ii(^)  =  f  u(x)e“'dx. 
Ja 


Using  Fourier  modes,  the  norm  on  the  space  L^{U)  can  be  characterized  by 

11/2 


'  r  _  ■ 

1/2 

f  o 

/  |u|  dx 

= 

/  |u|  di 

Ju 

Denote  the  Sobolev  space:  H'(f2)  =  {u  :  D^u  e  L^(Cl)  for  all  |s|  si  1}.  To  describe  the  boundary  operator  and  transparent 
boundary  condition  in  the  formulation  of  the  boundary  value  problem,  we  define  the  trace  functional  space 


|u  e  L^(U) :  y (1  +  <  ooj, 


whose  norm  is  defined  by 


IIh®(r)  ■ 


r 

(1 + . 

Jb 


by 


It  is  clear  that  the  dual  space  associated  with  H^(R)  is  the  space  H  ''(R)  with  respect  to  the  scalar  product  in  L^(R)  defined 


(u,  v)^  iivd^. 
Jb 


(2.7) 


To  simply  proofs,  we  shall  employ  positive  constants  C  and  C,  as  generic  constants  whose  precise  values  are  not  required 
and  may  be  changed  line  by  line  but  should  be  always  clear  from  the  context. 
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2.2.  Transparent  boundary  condition 

In  this  section,  we  shall  derive  a  boundary  operator,  which  maps  the  electric  field  to  its  normal  derivative,  and  introduce  a 
transparent  boundaiy  condition  on  the  aperture  of  the  cavity,  under  which  the  scattering  problem  may  be  reduced  into  a 
bounded  domain,  i.e.,  inside  the  cavity  fi,  from  the  unbounded  open  domain  fi  u  R^. 

By  taking  the  Fourier  transform  of  (2.5)  with  respect  to  x,  we  have  an  ordinary  differential  equation  with  respect  to  y: 

+  y>0.  (2.8) 

Since  the  solution  of  (2.8)  satisfies  the  radiation  condition  (2.6),  we  deduce  that  the  solution  of  (2.8)  has  the  analytical 
form 

=  u*(((,0)e‘"®^,  (2.9) 

f  (Kg  -  {2)1/2  foj.  1^1  < 

[  i({^  -  K-g)’/^  for  |{|  >  Ko. 

Taking  the  inverse  Fourier  transform  of  (2.9),  we  find  that 

u®(x,y)  =  f  u*({,0)e‘'’®-*'e“‘^i‘d{  in  R^. 

Jk 

Taking  the  normal  derivative  on  r‘^  u  F,  which  is  the  partial  derivative  with  respect  toy  on  F*^  u  F,  and  evaluating  aty  =  0 
yield 

anU*(x,y)|,=o=  [  (2.10) 

Ju 

where  n  is  the  unit  outward  normal  on  F"^  u  F,  i.e.,  n  =  (0, 1)^. 

For  any  given  u  on  F'^  u  F,  define  the  boundary  operator  T; 

Tu=  [  (2.11) 

Ja 

which  leads  to  a  transparent  boundary  condition  for  the  scattered  field  on  F*^  u  F: 

a„(u  -  u^if)  =  r(u  - 

Equivalently  it  can  be  written  as  a  transparent  boundary  condition  for  the  total  field 

dnU  =  Tu  +  g  onF'^uF,  (2.12) 

where 

g(x)  =  —  -2iKo  cos 

The  following  two  lemmas  are  concerned  with  the  continuity  and  analyticity  of  the  boundary  operator,  and  will  play  an 
important  role  in  the  proof  of  the  uniqueness  and  existence  of  the  solution  for  the  cavity  scattering  problem. 

Lemma  2.1.  The  boundary  operator  T  ■.  ^  continuous. 


where 

m)- 


Proof.  For  any  u,ve  it  follows  from  the  definitions  (2.7)  and  (2.11)  that 

(Fu,  v)  =  i  f  piiMi  =  i  f  /?({)(!  +  X  (1  +  {2)’/^u  x  (1  +  {2)’/''M{. 

Ju  Jr 


To  prove  the  lemma,  it  is  required  to  estimate 

mm 


—  00  <  {  <  00. 


Explicitly,  we  have 

\[({2_Kg)/(l+{2)]V2  for|{| 

It  can  be  verified  that  the  even  function  F({) 
culation  yields 


<  Ko, 

>  Kq. 

decreases  for  0  <  {  <  Ko  and  increases  for  Ko  <  {  <  oo.  Hence,  a  simple  cal- 
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F{i)  ^  max{F(0),  f(oo)}  =  max{K:o,  1}. 

Combining  above  estimates  and  using  the  Cauchy-Schwarz  inequality  yield 

l(ru,  v)\  ^  c||u||jji/2j[jj II 

where 

C  =  maxjfCo,  1}. 

Thus  we  have 

F“IIh->/^(r)  <  sup  C||u||„,/2,„).  □ 


Lemma  2.2.  Let  u  e  It  holds  that  Re(Tu,  u)  ^  0  and  Im(Tu,  u)  >  0.  Furthermore,  if  it  is  an  analytical  function  with  respect 

to  Re(Tu,  u)  =0  or  Im{Tu,  u)  =  0  implies  u  =  0. 


Proof.  By  definitions  (2.7)  and  (2.11),  we  find 

{Tu,u)  =  i  [  mfds. 

Ju 

Taking  the  real  part  gives 

Re{Tu.u)  =  -  /  ({^  -  fcg)’^^|u|^d{  s:  0 

7|{|>Ko 


and  taking  the  imaginary  part  yields 

lm(Tu,u)=  f  [kI  -  £^y^^\u\^di  ^  0. 

J\il<Ko 

Furthermore,  Re{Tu,  u)  =  0  implies  it  =  0  for  \^\  >  Ko  and  lm(Tu,  u)  =  0  implies  u  =  0  for  |{|  <  Kq.  If  it  is  assumed  to  be  an  ana¬ 
lytical  function  with  respect  to  c,  either  Re{Tu,u)  =  0  or  Im(Tu,u)  =  0  implies  u  =  0  for  e  R.  Taking  the  inverse  Fourier 
transform  of  u  =  0  yields  u  =  0,  which  completes  the  proof.  □ 


The  transparent  boundary  condition  is  derived  for  u  e  H'^^(R)  and  u  is  defined  over  u  F.  To  derive  a  transparent  bound¬ 
ary  condition  for  the  total  field  only  on  the  aperture  F,  we  need  to  make  the  zero  extension  as  follows;  for  any  given  u  on  F, 
define 


u(x) 


u  for  X  6  F, 

0  for  X  e  F". 


The  zero  extension  is  consistent  with  the  problem  geometry  where  the  cavity  is  placed  on  a  perfectly  conducting  ground 
plane  F"^,  i.e.,  the  total  field  u  is  required  to  be  zero  on  F‘^.  Based  on  the  extension  and  the  transparent  boundary  condition 
(2.12),  we  have  the  transparent  boundary  condition  for  the  total  field  on  the  aperture 

d„u  =  Tu+g  on  F.  (2.13) 


2.3.  Well-posedness 

We  now  present  a  variational  formulation  for  the  single  cavity  scattering  problem  and  give  a  proof  of  the  well-posedness 
of  the  boundary  value  problem. 

Define  a  trace  functional  space 

^{u-.ue  H’/^(R)}, 

whose  norm  is  defined  as  the  H’^^(R)-norm  for  its  extension,  i.e., 

ll“llHi/2(r)  ~  II“IIh'/2(r)'  (2.14) 

Define  a  dual  paring  (■,  ■)p  by 
(U,  v)y  =  J  uv. 

Obviously,  this  dual  paring  for  u  and  p  is  equivalent  to  the  scalar  product  in  L^(R)  for  their  extensions,  i.e., 

(U,  (U,  V). 
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by 


Denote  by  H  the  dual  space  of  i.e.,  H  ^^^(F)  ^  {H^^^(F))'.  The  norm  on  the  space  H  ^^^(F)  is  characterized 


llH-’/2(r)  ^ 


sup 


(U,  V}r 


sup 


(U,  V) 


eHV2(r)  II  ^llHV2(r)  reH'/2(R)  II  ^IIh’/2(r) 

Introduce  the  following  space 

H’(Q)  =  {u  6  H\a.)  :  u  =  0  on  S,  u|r  e 

which  is  clearly  a  Hilbert  space  with  the  usual  H'{Q)-norm. 

The  following  trace  regularity  in  is  useful  in  subsequent  analysis. 

Lemma  2.3.  For  any  u  e  it  holds  the  estimate 

ll“ll]^i/2(r)  ^  ‘^II“IIh'(£2)- 

where  C  is  a  positive  constant. 


(2.15) 


Proof.  Choose  a  positive  number  b  such  that  the  domain 

D  =  {(x,y)  E  :  -00  <x<oo,  -ti<y<0}  =  Rx  (-fa,  0) 

contains  the  cavity  fi,  i.e.,  fi  cc  D. 

Simple  calculation  yields 

fa|C(0)P  =  y^JC(y)|"dy  +  j ^|C(t)|^dfdy  <  J Ji:{y)fdy  +  b  J^^2\t:{y)\\!:'(y)\dy, 

which  implies  by  the  Cauchy- Schwarz  inequality  that 

(l  +  ia')’/'|C(0)P^C'(l  +  |iP)  fmfdy+  f  \C{y)fdy, 

J-b  J-b 

where  C  =  (1  +  fa“')'^^. 

Given  u  in  Hl(Q),  consider  the  zero  extension  to  the  domain  D: 

f  u  for  X  E  n, 
u(^)  =  t 

'  ^  lO  forxED\n. 

It  follows  from  the  definition  (2.14)  that 

l|u|l~./2(r)  =  II“IIh>/^(r)  =  ^(1  + 


Using  (2.16)  we  obtain 


(1  +  ^^)’/'|d({,0)|^  s;  C=(l  +  {2)  1°  |fi({,y)|'dy  + 


Iu'(i,y)l^dy4 


[(l+f")|fi({,y)|"  +  |u'(i,y)|"]dy. 


Noting 

II“IIh'{D)  [(''  +  {")l“(^-3')l"  +  |u'(^,y)|"]dyd^. 

Combining  (2.14),  (2.17)  and  (2.18)  yields 

ll“llHl/2(r)  ^  ^II“IIh'(D)  ^  *-II“IIh'(£!)i 

which  completes  the  proof.  □ 


(2.16) 


(2.17) 


(2.18) 


Multiplying  the  complex  conjugate  of  a  test  function  v  e  HJ(Q)  on  both  sides  of  (2.1 ),  integrating  over  f2,  using  the  inte¬ 
gration  by  parts  and  boundary  conditions  (2.2)  and  (2.13),  we  deduce  the  variational  formulation  for  the  single  cavity  scat¬ 
tering  problem:  find  u  s  Hj(n)  such  that 

aj(u,  v)  =  (g,  v)y  for  all  v  e  hI{Q.),  (2.19) 

where  the  sesquilinear  form  is 

a](u,v)^  /  (Vu  •  Vfa  -  K^ufa)  -  (Tu,  Z2)f,. 
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Noting  the  definition  for  the  dual  paring  and  using  the  scalar  product  in  L^(IR),  we  have  an  equivalent  sesquilinear  form 

(2.20) 


aj(u,p)^  I  (Vu-Vv-K^uv)  -  (Tu,v). 
Ja 


Theorem  2.1.  The  variational  problem  (2.19)  has  at  most  one  solution. 

Proof.  It  suffices  to  show  that  u  =  0  in  Q  if  g  =  0.  If  u  satisfies  the  homogeneous  variational  problem  (2.19),  then  we  have 

\\Vu\^-K^\uf-)-{Tu.u)  =  0. 


/( 


Taking  the  imaginary  part  of  the  above  identity  yields 


/ 


Imic^lul^  +  lm(rfi,  u)  =  0. 


Recall  =  co^Bfig,  Ime  >  0,  and  Im(ru.u)  >  0  from  Lemma  2.2,  we  get 
lm(Tu,u)  =  0. 

Since  u  has  a  compact  support  on  the  x-axis,  u  is  analytical  with  respect  to  (.  It  follows  from  Lemma  2.2  again  that  u  =  0.  The 
transparent  boundary  condition  (2.12)  yields  that  9„u  =  0  on  H  u  F.  An  application  of  Holmgren  uniqueness  theorem  yields 
u  =  0  in  A  unique  continuation  result  in  [21  ]  concludes  that  u  =  0  in  fi.  □ 

Theorem  2.2.  The  variational  problem  (2.19)  has  a  unique  weak  solution  in  Hl(Q)  and  the  solution  satisfies  the  estimate 
where  C  is  a  positive  constant. 

Proof.  Decompose  the  sesquilinear  form  (2.20)  into  Qi  =  On  -  au,  where 

aii(u,  ti)=  /  \7u  -Vv  -  {Tii,v) 

Ja 

and 

au(u,v)^  /  K^uv. 

Ja 

We  conclude  from  Lemma  2.2  and  Poincare  inequality  that  O]  is  coercive  from 

Rean(u,u)=  f  |Vu|^  -  Re(ru,  u)  >  f  |Vup  for  all  u  e  Hj  (LI) . 

Ja  Ja 

Next  we  prove  the  compactness  of  au-  Define  an  operator  Ri  :  L^{Q)  H^(Q)  by 

aii(/<'iu,  p)  =  au(u,  v)  for  all  v  e 
which  explicitly  gives 

[  VKjU  ■  Vv  -  (TKjU,  v)  ^  [  K^uv  for  all  e  H5(Q). 

Ja  Ja 

Using  the  coercivity  of  On  and  the  Lax-Milgram  Lemma,  it  follows  that 

||lLiU||fji(jj)  ^  C||u||j2jjjj.  (2.21) 

Thus  Ki  is  bounded  from  L^(f2)  to  H'(fi)  and  H'(f2)  is  compactly  imbedded  into  Hence  Ri  :  L^(Ci)  L^(Q)  is  a  com¬ 
pact  operator. 

Define  a  function  w  e  L^(f2)  by  requiring  w  e  Hj(Q)  and  satisfying 
au(w,  v)  ^  (g,  for  all  z' e  Hj)^). 

It  follows  from  Lemma  2.3  and  the  Lax-Milgram  Lemma  that 

II^IIh'{£!)  ^  C||g||j^-i/2jpj.  (2.22) 

Using  the  operator  Ri,  we  can  see  that  the  variational  problem  (2.19)  is  equivalent  to  find  u  e  such  that 
(I  -  I<j)u  =  w. 


(2.23) 
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It  follows  from  the  uniqueness  result  of  Theorem  2.1  and  the  Fredholm  alternative  that  the  operator  /  -  has  a  bounded 
inverse.  We  then  have  the  estimate  that 

^  C||w||i2(„).  (2.24) 

Combining  (2.21)-(2.24),  we  deduce  that 

II^IIh'(Q)  ^  ll^^l^llH'(n)  +  II'^IIh'(£!)  ^  C|I*^IIl^(£J)  +  II'^IIh'(!2)  ^  CII^IIhVO)  ^  C||g||jj-i/2jpj, 

which  completes  the  proof.  □ 


3.  Two  cavity  scattering 

To  address  the  general  multiple  cavity  scattering  problem,  we  begin  with  the  discussion  on  the  two  cavity  scattering 
problem.  The  two  cavity  scattering  problem  shares  the  same  features  with  the  general  multiple  cavity  scattering  problem, 
but  it  is  easier  to  present  the  major  ideas  in  the  proof  of  the  well-posedness  of  the  solution  for  the  multiple  cavity  scattering 
problem. 

3.1.  A  model  problem 

As  shown  in  Fig.  2,  two  open  cavity  fii  and  f22.  enclosed  by  the  aperture  F i  and  r2  and  the  walls  S]  and  S2,  are  placed  on  a 
perfectly  conducting  ground  plane  Above  the  flat  surface  {3/  =  0}  =  Fi  u  r2  U  F"^,  the  medium  is  assumed  to  be  homoge¬ 
neous  with  a  positive  dielectric  permittivity  Sq.  The  medium  inside  the  cavity  f2]  and  £^2  is  inhomogeneous  with  a  variable 
dielectric  permittivity  £i{x,y)  and  £2{x,y),  respectively.  Assume  further  that  Ej(x,y)  e  L°°(C1),  Refij  >  0,  Imej  >  0  for  j  =1,2. 

We  consider  the  TM  polarization,  where  the  time-harmonic  Maxwell  equations  are  reduced  to  the  two-dimensional 
Helmholtz  equation 

Au  +  k^u  =  0,  infliUfl2UR^,  (3.1) 

together  with  the  homogeneous  Dirichlet  boundary  condition 

u  =  0,  onF'^uSiUS2.  (3.2) 

Let  the  plane  wave  u‘  be  incident  on  the  cavities  from  above.  Due  to  the  interaction  between  the  incident  wave  and  the 
ground  plane  and  the  cavities,  the  total  field  u  is  assumed  to  be  consisted  of  the  incident  field  u‘,  the  reflected  field  uf  and  the 
scattered  field  uf  The  scattered  field  is  required  to  satisfy  the  radiation  condition  (2.6). 

To  reduce  the  scattering  problem  from  the  open  domain  £li  u  £22  U  into  the  bounded  domains  £2i  and  £22,  we  need  to 
derive  transparent  boundary  conditions  on  the  aperture  Fi  and  F2.  Rewrite  (3.1)  and  (3.2)  into  two  single  cavity  scattering 


problem: 

Aui  +  k:^Ui  =  0  in  £2],  (3.3) 

Ui  =  0  on  Si  (3.4) 

and 

Au2  +  Kju2^Q  in  £22,  (3.5) 

U2  =  0  on  S2,  (3.6) 


where  kJ  =  co^Sjiig,  j  =  1,2.  Clearly,  if  u  is  the  solution  of  (3.1)  and  (3.2),  and  Ui  and  U2  are  solutions  of  (3.3)  and  (3.4)  and 
(3.5)  and  (3.6),  respectively,  then  we  have 

iii  =  u|n,  and  U2  =  u|n^. 


T'  Ti  r  Ta  r 


5i 


^2 


52 


Fig.  2.  The  problem  geometry  for  the  two  cavity  scattering  problem.  Two  open  cavities  Hi  and  Q2.  enclosed  by  the  apertures  Fi  and  r2  and  the  walls  Si  and 
S2.  are  placed  on  a  perfectly  conducting  ground  plane  F*^. 
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Due  to  the  homogeneous  medium  in  the  upper  half  space  and  the  radiation  condition  (2.6),  the  scattered  field  still 
satisfies  the  same  ordinaiy  differential  equation  (2.9)  after  taking  the  Fourier  transform  with  respect  to  x,  and  thus  the  total 
field  u  satisfies  the  transparent  boundary  condition  (2.12),  which  is  now  written  as 

dnU  =  Tu  +  g  onr‘^uriUr2.  (3.7) 


For  Uj(x,0),  j 
Uj(x,0)  = 


=  1,2,  define  their  extensions  to  the  whole  x-axis  by 
f  Uj(x,0)  for  X  6  Fj, 

\  0  for  X  6  R  \  Fj. 


For  the  total  field  u(x,0),  define  its  extension  to  the  whole  x-axis  by 

fui  (x,  0)  for  X  e  Fi , 

U2(x,0)  forx6F2, 

0  for  X  e  FF 


It  follows  from  the  definitions  of  these  extensions  that  we  have 


ii  =  U]  +  U2  onF‘^uFiUF2. 

The  transparent  boundary  condition  (3.7)  can  be  written  as 


dnii  =  Tu+g  onF'^uFiUF2, 

(3.8) 

which  leads  to  the  transparent  boundary  conditions  for  Ui  and  U2: 

d„Uj  =  Tuj +Tu2 +g  onFi 

(3.9) 

and 

dnU2  =  TU2  +  Tu^+g  onF2. 

(3.10) 

As  we  can  see  from  (3.9)  and  (3.10),  the  boundary  conditions  for  Ui  and  U2  are  coupled  with  each  other,  which  is  the  major 
difference  between  the  single  cavity  scattering  problem  and  the  multiple  cavity  scattering  problem. 

The  following  lemma  is  analogous  to  Lemma  2.2  and  plays  an  important  role  in  the  analysis  of  the  uniqueness  and  exis¬ 
tence  for  the  solution  of  the  two  cavity  scattering  problem. 


Lemma  3.1.  Let  u,ve  H’^^(R).  It  holds  that 

Re((Tu,  u)  +  (T V,  v)  +  (Tu,  zi)  +  (Tzi,  u))  ^  0 

and 

lm((Tu,u)  +  (Tv,  v)  +  (Tu,  v)  +  {Tv,u))  >  0. 

Furthermore,  if  it  and  v  are  analytical  functions  with  respect  to  c,  either 
Re((Tu,  u)  +  (T V,  v)  +  (Tu,  v)  +  (T v,  u))  =  0 


or 

lm((Tu,  u)  +  (T V,  v)  +  (Tu,  v)  +  (Tv,u})  ^  0 

implies 

u  +  v  =  0. 


Proof.  By  definitions  (2.7)  and  (2.11),  we  have 

(Tu,u)  +  (Tv,v)  +  (Tu,v)  +  (Tv,u)  =  i  f  Il\\u\^  +  \v\^  +  2Re(uv)\d^  =  i  f  p\u+v\^di. 

Ja  i  i  Ju 

Taking  the  real  part  of  above  identity  gives 

Re{(Tu,u)  +  (Tv,  v)  +  (Tu,  v)  +  (Tv,u))  =  -  /"  (c^  -  v\^d^  ^  0 

7|{|>k„ 

and  taking  the  imaginary  part  of  above  identity  yields 

lm{(Tu,u)  +  (Tv,v)  +  (Tu,v)  +  (Tv,u))  =  f  (kI  -  +  vfd^  ^  0. 

i|{|<K0 
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Furthermore, 

Re((Tu,  u)  +  {Jv,  V)  +  {Tu,  v)  +  {Tv,  u))  =  0  (3.11) 

implies 

u  +  v  =  0  for  If  I  >  Ko 

and 

Im((ru,  u)  +  {Tv,  V)  +  {Tu,  v)  +  {Tv,  u))  =  0  (3.12) 

implies 

u  +  v  =  0  forjfl  <  Ko. 

If  it  and  v  are  assumed  to  be  analytical  functions  with  respect  to  f,  then  either  (3.11)  or  (3.12)  implies  that  ii  +  v  =  0  for  all 
f  6  R.  Taking  the  inverse  Fourier  transform  of  u  +  t'  =  0  yields  u  +  v  =  0,  which  completes  the  proof.  □ 


3.2.  Well-posedness 


u  = 


We  now  present  a  variational  formulation  for  the  two  cavity  scattering  problem  and  give  a  proof  of  the  well-posedness  of 
solution  for  the  boundaiy  value  problem. 

Denote  f2  =  fii  un2,r  =  ri  ur2,  and  S  =  Si  US2.  Let 

■  Ui  infii, 

,  U2  in  122- 

Define  a  trace  functional  space 
H’/2(r)  =  H’/2(ri)  X 
whose  norm  is  characterized  by 


'HV2(r) 


Hi/2(ri) 


IIU2I 


Hi/2(r2)' 


Denote  H  ’^^(F)  =  H  ’^^(Fi)  x  H  '^^(F2),  which  is  the  dual  space  of  HF2(r).  The  norm  on  the  space  H  ’^^(F)  is  charac¬ 
terized  by 


V2(r)  - 


1*^1  llH-'/2(r,)  ■ 


\\U2\\ 


Introduce  the  following  space 
Hj(n)  =  H’,(Qi)xH4(a2), 
which  is  a  Hilbert  space  with  norm  characterized  by 


ll^lluhn)  —  ll*^illH'(fii)  ^  ll^2|lH’(n2)' 


Multiplying  the  complex  conjugate  of  test  function  V\  e  (f2])  on  both  sides  of  (3.3),  integrating  over  fl],  and  using  the 
integration  by  parts  and  the  boundary  conditions  (3.4)  and  (3.9),  we  have 


/ 


(Vui  •  VZ2,  -  KjUiZ^i)  -  {Till,  Vi)  -  {Tii2,  Vi)  =  {g,  Vi}y^. 

Similarly,  regarding  (3.5),  (3.6),  and  (3.10),  we  have 

(VU2  •  Vt>2  -  Kj'^2V2)  -  {Till,  V2)  -  {Tii2,  V2)  =  {g,  V2}r.  for  all  V2  e  H’  (02). 


/ 

«/  Qp 


(3.13) 


(3.14) 


Adding  (3.13)  and  (3.14),  we  deduce  the  variational  formulation  for  the  two  cavity  scattering  problem:  find  u  e  IHIs(f2), 
where  Ui  =  u|ji^  e  Hs,  (f2])  and  U2  =  u\^^  e  Hs2(Ci2),  such  that  for  all  where  f]  =  e  Hs, (fli)  and  V2  = 

v\a2  6  Hs2{Ti2),  it  holds 

a2(u,  V)  =  {g,  Vi)j.^  +  {g,  V2)r2:  (3.15) 

with  the  sesquilinear  form 

a2{u,v)^[  {Vui-Vvi-KjuiVi)+  [  {Vu2-Vv2-kIu2V2)  -  {Tiii,Vi)  -  {Tu2,Vi)  -  {Tui,V2)  -  {Tu2,V2).  (3.16) 
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Theorem  3.1.  The  variational  problem  (3.15)  has  at  most  one  solution. 


Proof.  It  suffices  to  show  that  Ui  =  0  in  fii  and  U2  =  0  in  f22  if  g  =  0.  If  Ui  and  Uz  satisfy  the  homogeneous  variational  prob¬ 
lem  (3.15)  in  Q]  and  f22.  then  we  have 

[  (IVUil^  -  K]\U^\^)  +  [  (|VU2|'  -  Kl\U2f)  -  (TU:,U,)  -  (TU2,U,)  -  {TUuU2)  -  (TU2M2}  =  0. 

Ja,  Jq2 

Recall  Im8i  >  0  and  Ime2  3;  0.  Taking  the  imaginary  part  of  the  above  identity  yields 
Im((Tui,Ui)  +  (Tu2,Ui)  +  (Tui,U2>  +  {Tu2,U2))  =  0. 

Since  U]  and  U2  have  compact  supports  on  the  x-axis,  U]  and  U2  are  analytical  with  respect  to  Hence  we  have  from 
Lemma  3.1  that 

Ui  +  U2  =  h  =  0. 

The  transparent  boundary  condition  (3.8)  yields  that  9„u  =  0  on  P  u  Fi  u  r2.  An  application  of  Holmgren  uniqueness  the¬ 
orem  yields  u  =  0  in  A  unique  continuation  result  in  [21]  concludes  that  Ui  =  0  in  Hi  and  U2  =  0  in  H2-  □ 


Theorem  3.2.  The  variational  problem  (3.15)  has  a  unique  weak  solution  u  in  Hj]!!)  and  the  solution  satisfies  the  estimate 
ll^lliHVn)  ^  C||g||n-i/2(p), 
where  C  is  a  positive  constant. 


Proof.  The  proof  is  analogous  to  that  for  Theorem  2.2.  We  sketch  it  here.  Decompose  the  sesquilinear  form  (3.16)  into 
02  =  021  -  022,  where 


02(11,1')=  /  VUi-Vt'l+/  VU2  •  Vt'2  -  (FUi,  Vi)  -  (TU2,  t’l)  -  (Tui,  t’2)  -  (ru2,  %) 

Ja,  Ja2 

and 

a2{u,v)=  /  KiUiZii+  /  k\u2V2- 

We  conclude  from  Lemma  3.1  and  Poincare  inequality  that  021  is  coercive  from 

Rea2i(u,u)=  f  |Viii|^+  [  |Vu2|^  -  Re((Tu,,Ui)  +  (Tuj.Ui)  +  (Tu,,U2>  +  (ru2,U2)  >  f  |Viii|^+  f  |Vu2|^ 
J  CI2  t/ Q2 


>  C,||u,||^.(n^)  +  C2||U2|lHi{n,,  5=  C||u||^„,(n,  for  all  u  e  H’(n). 

Next  we  prove  the  compactness  of  022.  Define  an  operator  K2  :  L^(H)  ^  H^H)  by 
02i(/<'2U,  v)  =  022(0,  V)  for  all  v  e  IHlJ(H). 

The  continuity  of  the  sesquilinear  form  O22  follows  from  the  Cauchy-Schwarz  inequality 

1022(0,0)1  <  Cl  ||Ol  Ol  +  C2||02||j2(jj^j||  02||^2jjj^j  ^  C||u||,^2(jj)||0||^2(J2), 

where  L^(H)  =  L^{Hi)  x  L^{Cl2)-  Using  the  Lax-Milgram  Lemma  and  the  continuity  of  O22,  we  obtain 

||i^20|||nq;2)  ^  C||o|||^2(;2).  (3-17) 

Thus  K2  is  bounded  from  L^(H)  to  IHl'(H)  and  IHI^(H)  is  compactly  imbedded  into  L^(H).  Hence  R2  :  L^{H)  ^  L^{H)  is  a  com¬ 
pact  operator. 

Define  a  function  w  e  L^(H)  by  requiring  w  e  IHl|(f2)  and  satisfying 
a2i(w,  v)  =  (g,  Oi)p,  +  (g,  02>r,  for  all  v  e  H’(n). 

It  follows  from  Lemma  2.3  and  the  Lax-Milgram  Lemma  that 


l|w|lH'(n)  ^  C||g||„-i/2(p).  (3.18) 

Using  the  operator  K2.  we  can  see  that  the  variational  problem  (3.15)  is  equivalent  to  find  u  s  L^{H)  such  that 

{1-K2)u  =  w.  (3.19) 


It  follows  from  the  uniqueness  result  of  Theorem  3.1  and  the  Fredholm  alternative  that  the  operator  /  -  K2  has  a  bounded 
inverse.  We  then  have  the  estimate  that 
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<  C||w||L2(n,. 

From  (3.17)-(3.20),  we  deduce  that 

II^ILvn)  ^  II^^2K|Ih'(!J)  +  II^IL'(n)  ^  C||u||,^2(J2) 

which  completes  the  proof.  □ 


II'^IIh'{£!)  ^  ^II'^IIh'(£!)  ^  C|IS|Ih-'/2(p), 


(3.20) 


4.  Multiple  cavity  scattering 

Now  we  generalize  the  model  problem  and  techniques  for  the  two  cavity  scattering  to  the  case  of  multiple  cavity  scatter¬ 
ing.  The  proofs  are  analogous  to  those  for  the  two  cavity  scattering  problem.  For  completeness,  we  shall  briefly  discuss  the 
model  problem  and  the  solution  for  the  multiple  scattering  problem. 

4.1.  A  model  problem 


As  shown  in  Fig.  3,  we  consider  a  situation  with  n  cavities,  where  the  multiple  open  cavity  fli, . . .  ,f2„,  enclosed  by  the 

aperture  Fi, _ r„  and  the  walls  Si _ ,S„,  are  placed  on  a  perfectly  conducting  ground  plane  F"^.  Above  the  flat  surface 

{y  =  0}  =  F,  U  ■  ■  ■  U  F„  u  F'^,  the  medium  is  assumed  to  be  homogeneous  with  a  positive  dielectric  permittivity  Eq.  The  med¬ 
ium  inside  the  cavity  Qj  is  inhomogeneous  with  a  variable  dielectric  permittivity  £;(x,y).  Assume  further  that 
6  Refij  >  0,  ImSj  ^  0  for  j  =  1, . . . ,  n. 

We  consider  the  same  model  of  the  two-dimensional  Helmholtz  equation  for  the  total  held: 

Au  +  k^u  =  0,  in  u  •  •  •  U  U  R^,  (4.1) 

together  with  the  perfectly  electric  conductor  condition 

u  =  0,  on  F'^  U  Fi  U  •  •  ■  U  F„.  (4.2) 

The  total  held  u  is  assumed  to  consist  of  the  incident  field  u‘,  the  reflected  field  uf  and  the  scattered  held  uf  where  the 
scattered  field  is  required  to  satisfy  the  radiation  condition  (2.6). 

To  derive  the  transparent  boundary  condition  on  the  aperture  Fj,  j  =  1 , . . . ,  n,  we  reformulate  the  multiple  cavity  scatter¬ 
ing  problem  (4.1)  and  (4.2)  into  n  single  cavity  scattering  problems  which  are  coupled  through  the  boundaiy  condition. 

AUj  +  Kfuj  =  0  in  Qj,  (4.3) 


Uj  =  0  on  Fj, 


where  k?  =  ocAEjji^,  j  =  1,. . .  ,n. 

For  Uj(x.O).  j  =  1 . n,  define  its  zero  extension  to  the  whole  x-axis. 


Uj(x,0)  = 


u{x,  0)  = 


'  Uj{x,  0)  for  X  e  Fj, 

0  for  X  6  R  \  Fj. 

For  the  total  field  u(x.O),  define  its  extension  to  the  whole  x-axis  by 
Uj(x,0)  forxeFj, 

0  forxEFF 

It  follows  from  the  definition  of  the  extension  that  we  have 


u  =  '^Uj  on  F‘^  u  Fi  u  •  •  •  U  F„. 

j=i 


(4.4) 


pC 

Ti 

pC 

pC  pC 

r„ 

pC 

Hi 

^2 

^2  S„ 


Fig.  3.  The  problem  geometry  for  the  multiple  cavity  scattering  problem.  Multiple  open  cavities  Hi , . . . ,  Hn .  enclosed  by  the  apertures  Fi , . . . ,  Fn  and  the 
walls  Si , . . . ,  Sn,  are  placed  on  a  perfectly  conducting  ground  plane  F^ 
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The  transparent  boundary  condition  can  be  written  as 
dnU  =  Tu+g  on  r'^  U  Fi  U  •  •  •  U  r„, 
which  leads  to  the  transparent  boundary  condition  for  uf. 


d„Uj  =  Tuj  +  ^TUi  +g  on  Fj. 


Lemma  4.1.  Let  Uj  e  j  =  1 _ ,n.  It  holds  that 

Re^^{TUj,u,)  ^  0 

i=i  j=i 


and 


5=  0. 

i=i  j=i 

Furthermore,  if  iij.  j  =  I _ _  n  are  analytical  functions  with  respect  to  either 

ReX:i(Tu;,u.)  =  0 

i=l  ;=1 


1=1  j=i 


implies 


=  0- 
j=i 


Proof.  By  definitions  (2.7)  and  (2.11),  we  have 


^^(Fuj,u,.)  =  i  [  f(^) 

j=l  j=l 


Taking  the  real  part  gives 


j=i 


di 


Re^^(TUj,  u,)  =  -  (  -  kI) 

1-1  j=:l 

and  taking  the  imaginary  part  yields 

Im^^(ruj,u,)  =  /  (Kg-c") 

i=l  j=l  '^kK^'o 


1=1 


d^^O 


Yuj 
1=1 


^  0. 


Furthermore, 

n  n 


Re^^(Tu;,u,)  =  0 

i=l  ;=1 


implies 


Ydj  =  0  for  1^1  >  Ko, 
1=1 


and 


Im^^(FUj,u,)  =  0 

i=i  j=i 


(4.5) 


(4.6) 


(4.7) 


(4.8) 
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implies 

=  0  for  1^1  <  Kg. 
j=i 

If  Uj,  j  =  1, . . .  .n,  are  assumed  to  be  analytical  functions  with  respect  to  i,  then  either  (4.7)  or  (4.8)  implies  that 

n 

y^Uj  =  0  for  all  {  6  R.  (4.9) 

j=i 

The  proof  is  completed  by  taking  the  inverse  Fourier  transform  of  (4.9).  □ 


4.2.  Well-posedness 


We  now  present  a  variational  formulation  for  the  multiple  cavity  scattering  problem  and  sketch  the  proof  for  the  well- 
posedness  of  the  boundary  value  problem. 

Denote  =  r  =  riU-urn,  and  S  =  Si  u  ■  ■  ■  U  S„.  For  simplicity,  we  shall  use  the  same  notation  as  those 

adopted  in  Section  3  for  the  two  cavity  scattering  problem.  The  following  notation  for  the  general  multiple  cavity  scattering 
problem  are  actually  consistent  with  those  for  the  two  cavity  scattering  problem  when  letting  n  =  2. 

Define  a  trace  functional  space 

H’/yr)  =  H’/yri)  X  •  •  •  X  H’/yr„). 

Its  norm  is  characterized  by 


'HV2(r) 


j=i 


HV2(ri)' 


Denote  H  ’/^(r)=H  ’/^(Fi)  x  ■  ■  ■  x  H  ’^^(r„),  which  is  the  dual  space  of  IH|F2(r).  The  norm  on  the  space  H  '/^(F)  is  char¬ 
acterized  by 

n 

j=i 

Introduce  the  following  space 

Hj(n)  =  H’yQi)x...xH’„(a„), 


which  is  a  Hilbert  space  with  norm  characterized  by 

n 

II^ILfin)  = 

j=i 

Multiplying  the  complex  conjugate  of  test  function  Vj  e  Hl(Qj)  on  both  sides  of  (4.3),  integrating  over  Qj,  and  using  the 
integration  by  parts  and  boundary  conditions  (4.4)  and  (4.6),  we  obtain 

f  {VUj  ■  VVj  -  KjujVj)  -  ^(Tu„  Vj)  =  (g,  (4.10) 

JQj 

Taking  summation  of  (4.10)  forj  ^  1, . . .  ,n,  we  deduce  the  variational  formulation  for  the  multiple  cavity  scattering  prob¬ 
lem:  find  u  e  HlJ(Q)  with  uj  ^  such  that 


n 

a3(u,  v)  =  y^(g,  Vj)Y^  for  all  v  a  HlJ(n), 
j=i 

where  the  sesquilinear  form 

a3(u,  =  E  /  “  EE<’^“'’ 

j=l  j=l  i=l 


(4.11) 


(4.12) 


Theorem  4.1.  The  variational  problem  (4.11)  has  at  most  one  solution. 
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Proof.  It  suffices  to  show  that  Uj  ^  0  in  for  j  ^  1, . . . ,  n  if  g  ^  0.  If  Uj  satisfy  the  homogeneous  variational  problem  in  Qj, 
then  we  have 


f)  [  (\\7uf  -  Kf\u/)  -  ^^(Tu„Uj)  =  0. 

j=l  j=\  {=1 


Noting  \mSj  ^  0  and  taking  the  imaginary  part  yields 

Im^^(rUi,Uj)  =  0. 

H  i=l 

Since  u,  has  a  compact  support  on  the  x-axis,  Uj  is  analytical  with  respect  to  Hence  we  have  from  Lemma  4.1  that 
u  =  =  0. 

j=i 


By  the  definition  of  the  extensions  iij  we  obtain 
Uj  ^0  on  r'^  U  Fi  U  •  •  •  U  r„. 

The  transparent  boundary  condition  (4.5)  yields  that  =  0  on  F'^  u  Fi  u  ■  ■  ■  U  Fn.  An  application  of  Holmgren  uniqueness 
theorem  yields  u  =  0  in  A  unique  continuation  result  in  [21  ]  concludes  that  Uj  =  0  in  fij  for  j  =  1 , . . . ,  n.  □ 


We  have  the  following  well-posedness  result  for  the  general  multiple  cavity  scattering  problem.  The  proof  is  similar  in 
nature  as  that  of  the  two  cavity  model  problem  and  is  omitted  here  for  brevity. 

Theorem  4.2.  The  variational  problem  (4.11)  has  a  unique  weak  solution  u  in  H|(f2)  and  the  solution  satisfies  the  estimate 

II^IL'(n)  ^  C||SIL-V2(r)> 

where  C  is  a  positive  constant. 


5.  Numerical  experiments 

In  this  section,  we  discuss  the  computational  aspects  for  solving  the  multiple  scattering  problems,  including  the  finite 
element  approximation,  a  block  Gauss-Seidel  iteration  method  for  the  discrete  weak  formulation,  and  an  alternative  trans¬ 
parent  boundary  condition. 

5.1.  Finite  element  formulation 

Let  Mj  be  a  regular  conforming  triangulation  of  Qj  with  Mj  small  triangular  element  and  Vj  c  be  the  conforming  linear 
finite  element  space  over  Mj.  Denote  V  =  Vi  x  . . .  x  F„.  The  finite  element  approximation  to  the  multiple  cavity  scattering 
problem  (4.11)  is:  find  u’'  eV  with  uj'  e  Vj  such  that 

n 

a3(u'',  v^)  =  for  all  e  V,  (5.1) 

j=i  ' 

where  the  sesquilinear  form 

n  n  n  n 

a3(u\  v'')  =  Y,  -  Kfu^vf)  -  uj').  (5.2) 

j=l  j=l  i=l 

For  any  1  ^  j n,  we  denote  by  Vj  the  set  of  N,  vertices  of  Mj  which  are  not  on  the  cavity  wall  Sj,  and  let  (Pj{r)  e  Vj  be  the 
nodal  basis  function  belonging  to  vertex  r  e  Vj.  Using  the  basis  functions,  the  solution  of  (5.1)  is  represented  as 

=  5]“j(r)(?>j(r)- 

reV,’ 

The  discrete  problem  (5.1)  is  equivalent  to  the  following  system  of  algebraic  equations 


Ai-Bi,i  -62,1 

— Bi_2  A2  —  B22 

-Bn.l 

—Bn.2 

u, 

U2 

■gr 

g2 

-Bin  -B2,n 

An  —  Bn,n . 

.Un. 

W  .. 

1 _ 

(5.3) 
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where  each  U/  is  an  unknown  vector  whose  entries  are  Uj{r)  ^  uj‘(r)  for  all  r  g  V/,  Aj  is  the  stiffness  matrix  for  the  discrete 
problem  of  i/|  and  its  entries  are  defined  by 

A(r,r')  =  ^  [V(pj(r)  ■  V(Pj{r')  -  Kj(pj(r)(p^ir')j  for  all  r,!^  e  Vj, 
the  entries  of  B,j  are  defined  by 


Bj, i(r,  r')  =  {T^i{r),  <pj(r'))  for  all  r,  r'  e  Vj  n  Tj 


and  the  entries  of  each  vector  gj  are  given  by 
g,-(r)  =  {g,  (pj{T))^^  for  all  r  e  Vj  n  r^. 

The  linear  system  (5.3)  requires  to  solve  a  large  and  coupled  equation,  particularly  for  multiple  cavities.  We  present  an 
efficient  block  Gauss-Seidel  method,  which  may  be  written  as  follows:  given  (uj”*,  ■  ■  ■  define  (uj*'*, . . .  >  1  by 

the  solution  of  the  following  system  of  equations 


j-i 


(A  -  %)uf = g; + +  E 


(H) 


1  s:  n. 


N+i 


(5.4) 


The  block  Gauss-Seidel  iteration  (5.4)  is  equivalent  to  apply  the  finite  element  method  for  solving  the  following  problem: 
Let  . . . ,  =  (0, . . . ,  0)^.  define  . . . ,  for  k  >  1  by  the  solutions  of  the  following  decoupled  equations  for 

j=l,...,n: 

Auf'>  +  KfuP  =  0  in  £ij,  (5.5) 

up  =  0  on  Sj,  (5.6) 

dnuf  =  +g  on  Fj,  (5.7) 

i=l  i=j+l 


Therefore,  we  only  need  to  solve  a  well-posed  single  cavity  scattering  problems  (5.5)-(5.7)  for  the  block  Gauss-Seidel 
method  at  each  iteration. 


5.2.  Transparent  boundary  condition 

Based  on  the  boundary  operator  (2.11),  the  transparent  boundary  conditions  (2.12),  (3.7),  and  (4.5)  are  useful  to  carry 
mathematical  analysis  for  the  boundary  value  problems.  However,  they  are  not  convenient  to  be  implemented  numerically. 
In  practice,  we  adopt  an  alternative  and  equivalent  transparent  boundary  condition  [38], 

Let 

G(x.x')  =  L  [h<’>(Ko/j)  -  hP(KoP)\ 

be  the  Green  function  of  the  Helmholtz  equation  in  the  upper  half  space,  where  is  the  Hankel  function  of  the  first  kind 
with  order  zero,  X  =  (x,y),x'  =  (x',y'),p  =  |x-  x'|,p  =  |x  -  x'|,  and  x'  =  (x',  -y')  is  the  image  ofx'  with  respect  to  the  real  axis. 
By  the  Green’s  theorem  and  the  radiation  condition,  we  obtain 

dyU^x,  0)  =  ^  ^  -  x'|)u^(x',  0)dx', 

where  Hi(l)  is  the  Hankel  function  of  the  first  kind  with  order  one.  Hence,  the  alternative  boundary  condition  is 

d„u  =  Tu  +  g  on  r,  (5.8) 

where  the  boundary  operator  T  is  defined  as 

ru  =  ^  ^  ^^HS”(Ko|x-x'|)u(x',0)dx'.  (5.9) 

Here  the  integral  is  understood  in  the  sense  of  Hadamard  finite-part.  For  multiple  cavities  with  apertures  given  as 
F]  u  . . .  u  r„,  the  boundary  operator  is  defined  as 

=  P^^H'’>(Ko|x-x'|)u(x',0)dx'. 
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The  boundary  operator  (5.9)  can  be  approximated  by 

m 

Tu(Xi)  fs  Y,SikU{Xk,0), 

k=l 

where 


and 


Kofl.J,(Ko|X,-Xfc|) 


tik  = 


S(l-'n2) 
_  2_ 
ht 


1  In  1'-''! 


for  \i-k\  =  1, 
for  |i  -  fe|  =  0, 
for  |i  -  k\  ^  2, 


where  h*  is  the  step  size  of  the  partition  for  the  cavity  aperture  F.  Therefore,  the  boundary  integral  (Tu,  v)  in  the  weak  for¬ 
mulation  for  the  cavity  scattering  problem  can  be  approximated  by  any  numerical  quadratures. 


5.3.  Numerical  examples 

In  this  section,  two  examples  of  the  cavity  scattering  model  problems  are  presented  to  illustrate  the  proposed  method. 
The  physical  parameter  of  interest  is  the  radar  cross  section  (RCS),  which  is  defined  by 

cr-^\P{(pf. 

Ko 

Here  iji  is  tbe  observation  angle  and  P  is  the  far-field  coefficient  given  by 

P((p)  =  ^  sin  (/)  J  u(x,  0)e‘'^'“'‘“®''’dx. 

When  the  incident  and  observation  directions  are  tbe  same,  a  is  called  the  backscatter  RCS,  which  is  defined  by 
Backscatter  RCS((p)  =  101og,Q(T((p)dB. 

In  the  following  two  examples,  the  numerical  results  are  obtained  by  using  a  linear  finite  element  over  triangles  at  the 
wavenumber  Kg  =  n,  which  accounts  for  the  wavelength  X  =  2. 

Example  1.  Consider  a  plane  wave  scattering  from  the  model  problem  of  three  identical  rectangular  cavities  at  normal 
incidence,  i.e.,  0  =  0.  Each  cavity  is  0.5/  wide  and  0.5/  deep;  there  is  0.5/  distance  away  from  each  other.  Specifically,  the 
domains  are  given  as  follows  for  the  three  rectangular  cavities: 

Cavity!:  [-2.5, -1.5]  x  [-1.0, 0.0], 

Cavity2:  [-0.5, 0.5]  x  [-1.0, 0.0], 

CavityS:  [1. 5, 2.5]  x  [-1.0, 0.0]. 

Each  rectangular  domain  is  divided  into  100  x  100  small  equal  rectangles  and  then  each  small  rectangle  is  subdivided  into 
two  equal  triangles. 

Table  1  lists  the  six  types  of  combination  of  cavities,  where  £  =  1  means  an  empty  or  unfilled  cavity  and  £  =  1  +  i  stands 
for  a  filled  cavity  with  some  homogeneous  absorbing  medium. 

Eigs.  4-9  show  the  magnitude  and  the  phase  of  the  total  field  on  the  apertures  at  the  normal  incidence  and  the  backscat¬ 
ter  RCS  for  all  the  six  types  of  cavities.  These  numerical  results  are  obtained  by  using  the  block  Gauss-Seidel  iterative  meth¬ 
od.  To  show  the  convergence  of  the  iterative  method,  we  define  the  error  between  two  consecutive  approximations 


Table  1 

Six  types  of  combination  of  cavities  for  Example  1. 


e 

Type  1 

Type  2 

Type  3 

Type  4 

Type  5 

Type  6 

Cavity  1 

1.0 

1.0 

1.0 

1.0 -ri 

1.0-ri 

1.0  +  i 

Cavity  2 

1.0 

1.0 

1.0  +  i 

1.0 -ri 

1.0 

1.0  +  i 

Cavity  3 

1.0 

1.0 -ri 

1.0 

1.0 

1.0  +  i 

1.0  +  i 
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X  X  Observation  angle  (degree) 

Fig.  4.  The  magnitude  and  the  phase  of  the  total  field  on  the  aperture  at  the  normal  incidence,  and  the  backscatter  RCS  for  Example  1  with  type  1  cavity: 
(left)  magnitude;  (middle)  phase;  (right)  backscatter  RCS. 


X  X  Observation  angle  (degree) 

Fig.  5.  The  magnitude  and  the  phase  of  the  total  field  on  the  aperture  at  the  normal  incidence,  and  the  backscatter  RCS  for  Example  1  with  type  2  cavity: 
(left)  magnitude;  (middle)  phase;  (right)  backscatter  RCS. 


Observation  angle  (degree) 


Fig.  6.  The  magnitude  and  the  phase  of  the  total  field  on  the  aperture  at  the  normal  incidence,  and  the  backscatter  RCS  for  Example  1  with  type  3  cavity: 
(left)  magnitude;  (middle)  phase;  (right)  backscatter  RCS. 


efc  =  max{||uf>  -  uj 


(H) 


where  k  is  the  number  of  iteration.  Fig.  10  shows  the  error  of  the  block  Gauss-Seidel  method  for  two  consecutive  approx¬ 
imations  again  the  number  of  iterations  for  all  three  types  of  cavities.  It  can  be  seen  from  Fig.  10  that  the  type  1  cavity  takes 
the  most  number  of  iterations  to  reach  the  same  level  of  accuracy  as  the  other  types  of  cavities;  the  type  6  uses  the  least 
number  of  iterations  to  reach  the  same  level  of  accuracy;  the  types  2  and  3  and  the  types  4  and  5  take  almost  the  same  num¬ 
ber  of  iterations  to  reach  the  same  level  of  accuracy,  respectively,  while  the  types  4  and  5  take  less  number  of  iterations  than 
the  types  2  and  3.  The  reason  is  that  the  types  2-6  cavities  are  filled  with  complex  medium  which  accounts  for  the  absorp¬ 
tion  of  the  energy  and  thus  the  damping  of  the  amplitude  of  the  field;  the  more  number  of  cavities  are  filled  with  the  absorp¬ 
tion  material,  the  less  number  of  iterations  are  needed  to  reach  certain  level  of  accuracy. 
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Fig.  7.  The  magnitude  and  the  phase  of  the  total  field  on  the  aperture  at  the  normal  incidence,  and  the  backscatter  RCS  for  Example  1  with  type  4  cavity: 
(left)  magnitude;  (middle)  phase;  (right)  backscatter  RCS. 
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Fig.  8.  The  magnitude  and  the  phase  of  the  total  field  on  the  aperture  at  the  normal  incidence,  and  the  backscatter  RCS  for  Example  1  with  type  5  cavity: 
(left)  magnitude;  (middle)  phase;  (right)  backscatter  RCS. 
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Fig.  9.  The  magnitude  and  the  phase  of  the  total  field  on  the  aperture  at  the  normal  incidence,  and  the  backscatter  RCS  for  Example  1  with  type  6  cavity: 
(left)  magnitude;  (middle)  phase;  (right)  backscatter  RCS. 


Example  2.  Consider  the  scattering  of  a  quintuple  cavity  model.  Let  a  plane  wave  be  incident  onto  five  identical  rectangular 
cavities  at  the  normal  direction.  Each  cavity  is  0.5/t  wide  and  0.5/  deep;  there  is  0.5/  distance  away  from  each  other.  The 
domains  for  the  five  rectangular  cavities  are  given  as  follows: 

Cavity!:  [-4.5, -3.5]  x  [-1.0, 0.0], 

Cavity  2:  [-2.5, -1.5]  x  [-1.0, 0.0], 

Cavity  3:  [-0.5,0.5|  x  [-1.0,0.0[, 

Cavity  4:  [1.5, 2.5]  x  [-1.0, 0.0], 

Cavity  5:  [3.5,4.5[  x  [-1.0,0.0[. 


P.  Li,  A.  Wood/Joumal  of  Computational  Physics  240  (2013)  100-120 


119 


Fig.  10.  The  error  of  two  consecutive  approximations  again  the  number  of  iterations  for  all  six  types  of  cavities  in  Example  1. 
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Fig.  11.  The  magnitude  and  the  phase  of  the  total  field  on  the  aperture  at  the  normal  incidence,  and  the  backscatter  RCS  for  Example  1:  (left)  magnitude; 
(middle)  phase;  (right)  backscatter  RCS. 


Again,  each  rectangular  domain  is  divided  into  100  x  100  small  equal  rectangles  and  then  each  small  rectangle  is  subdivided 
into  two  equal  triangles.  The  cavity  2  and  4  are  empty  cavities  with  e  =  1.0,  and  the  cavities  1,  3,  and  5  are  filled  with  the 
same  homogeneous  medium  with  e  =  1 .0  +  i.  Fig.  1 1  shows  the  magnitude  and  the  phase  of  the  total  field  on  the  apertures  at 
the  normal  incidence  and  the  backscatter  RCS. 


6.  Concluding  remarks 

The  problem  of  electromagnetic  scattering  by  cavities  embedded  in  the  infinite  two-dimensional  ground  plane  is  an 
important  area  of  research.  This  paper  marks  the  first  known  mathematical  investigation  into  the  phenomenon  of  scattering 
by  multiple  cavities.  We  solve  the  problem  by  reducing  the  overall  scattering  problem  to  coupled  single  cavity  scattering 
problems  via  the  introduction  of  a  novel  transparent  boundaiy  condition  over  the  cavity  apertures.  Uniqueness  and  exis¬ 
tence  of  the  variational  formulation  for  the  multiple  cavity  scattering  problem  is  achieved.  Numerical  experiments  of  model 
problems  demonstrate  the  efficiency  and  accuracy  of  our  numerical  methods.  Our  numerical  approach  gives  rise  to  the  wide 
spread  applicability  of  numerical  solvers  developed  for  single  cavity  models  to  the  general  multiple  cavity  setting.  Numerical 
examples  indicate  that  the  convergence  of  the  Gauss-Seidel  iterative  method  depends  on  the  wavenumber  Kq  and  the 
dielectric  permittivity  fij  inside  the  cavities;  it  does  not  depend  on  the  separation  distance  between  cavities  even  for  the  dis¬ 
tance  which  is  smaller  than  half  the  wavelength.  When  the  number  of  cavities  is  large,  it  does  take  a  lot  more  time  for  the 
method  to  go  over  all  the  cavities  in  order  to  account  for  the  multiple  interaction  among  them.  Future  research  includes  the 
rigorous  analysis  of  the  convergence  of  the  Gauss-Seidel  iterative  method,  which  requires  the  stability  analysis  of  the  cavity 
scattering  problem  (see  [12]  for  a  single  cavity  model).  Multiple  over-filled  cavity  models  are  another  natural  consideration. 
Finally,  we  note  that  transient  scattering  problems  involving  multiple  cavities  can  be  particularly  interesting  and 
challenging. 
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